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Abstract 

There are reasons to believe that mechanisms exist in the solar interior which lead 
to random density perturbations in the resonant region of the Large Mixing Angle 
solution to the solar neutrino problem. We find that, in the presence of these density 
perturbations, the best fit point in the {siv? 26 , A.m?') parameter space moves to 
smaller values, compared with the values obtained for the standard LMA solution. 
Combining solar data with KamLAND results, we find a new compatibility region, 
which we caU VERY-LOW LMA, where sin^ 26^ ^ 0.6 and Am^ 2 x 10"^ eV^, 
for random density fluctuations of order 5% < ^ < 8%. We argue that such values 
of density fluctuations are still allowed by helioseismological observations at small 
scales of order 10 - 1000 km deep inside the solar core. 
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Assuming CPT invariance, the electronic antineutrino z/g disappearance as well 
as the neutrino energy spectrum observed in KamLAND [1] are compatible 
with the predictions based on the Large Mixing Angle (LMA) realization of 
the MSW mechanism, resonantly enhanced oscillations in matter [2,3]. This 
compatibility makes LMA the best solution to the solar neutrino anomaly [4]- 
[10]. The best fit values of the relevant neutrino parameters which generate 
such a solution are Am? = 7.3 x 10~^ eV^ and tan^ 6 = 0.41 with a free 
boron neutrino flux = 1.05 [11] with the 1 a intervals: Am^ = (6.2 ^ 
8.4) X 10-5 eV^ and tan^^ = 0.33 ^ 0.54. 
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KamLAND results exclude not only other MSW solutions, such as the small 
mixing angle realization, but also the just-so solution and the alternative solu- 
tions [12] to the solar neutrino problem based on non-standard neutrino inter- 
actions [13]-[15], resonant spin-flavor precession in solar magnetic fleld [16]-[19] 
and violation of the equivalence principle [20]-[23]. Therefore, the LMA MSW 
solution appears, for the first time in the solar neutrino problem history, as 
the unique candidate to explain the anomaly. 

Such an agreement of the LMA MSW predictions with the solar neutrino 
data is achieved assuming the standard approximately exponentially decaying 
solar matter distribution [24]- [26]. This prediction to the matter distribution 
inside the sun is very robust since it is in good agreement with helioseismology 
observations [27]. 

When doing the fitting of the MSW predictions to the solar neutrino data, it 
is generally assumed that solar matter do not have any kind of perturbations. 
I.e., it is assumed that the matter density monotonically decays from the 
center to the surface of the Sun. There are reasons to believe, nevertheless, 
that the solar matter density fiuctuates around an equilibrium profile. Indeed, 
in the hydro-dynamical approximation, density perturbations can be induced 
by temperature fiuctuations due to convection of matter between layers with 
different temperatures. Considering a Boltzman distribution for the matter 
density, these density fluctuations are found to be around 5% [28]. Another 
estimation of the level of density perturbations in the solar interior can be 
given considering the continuity equation up to flrst order in density and 
velocity perturbations and the p-modes observations. This analysis leads to a 
value of density fluctuation around 0.3% [29]. 

The mechanism that might produce such density fluctuations can also be 
associated with modes excited by turbulent stress in the convective zone [30] 
or by a resonance between g-modes and magnetic Alfven waves [31]. As the 
g-modes occur within the solar radiative zone, these resonance creates spikes 
at speciflc radii within the Sun. It is not expected that these resonances alter 
the helioseismic analyses because as they occur deep inside the Sun, they do 
not affect substantially the observed p-modes. This resonance depends on the 
density proflle and on the solar magnetic fleld, and as mentioned in Ref. [31], 
for a magnetic fleld of order of 10 kG the spacing between the spikes is around 
100 km. In the analyses presented in Ref. [31] the values considered for the 
magnetic fleld are the ones that satisfy the Chandrasekar limit, which states 
that the magnetic fleld energy must be less than the gravitational binding 
energy. 

Considering helioseismology, there are constraints on the density fluctuations, 
but only those which vary over very long scales, much greater than 1000 
km [32,33,34]. In particular, the measured spectrum of helioseismic waves is 
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largely insensitive to the existence of density variations whose wavelength is 
short enough - on scales close to 100 km, deep inside within the solar core - to 
be of interest for neutrino oscillation, and the amplitude of these perturbations 
could be large as 10% [31]. 

So, there is no reason to exclude density perturbations at a few percent level 
and there are theoretical indications that they really exist. 

In the present paper, we study the effect on the Large Mixing Angle parameters 
when the density matter fluctuates around the equilibrium profile. We consider 
the case in which these fluctuations are given by a random noise added to an 
average value. This is a reasonable case, considering that in the lower frequency 
part of the Fourier spectrum, the p-modes resembles that of noise [28]. Also, 
considering the resonance of g-modes with Alfven waves, the superposition 
of several different modes results in a series of relatively sharp spikes in the 
radial density profile. The neutrino passing through these spikes fell them as a 
noisy perturbation whose correlation length is the spacing between the density 
spikes [31]. 

Noise fiuctuations have been considered in several cases. Loretti and Bal- 
antckin [35] have analyzed the effect of a noisy magnetic field and a noisy 
density. In particular, for a noisy density, they found that the MSW transi- 
tion is suppressed for a Small Mixing Angle. For the case of strongly adiabatic 
MSW transitions and large fluctuations, the averaged transition probability 
saturates at one-half. Refs. [28] and [36] have analyzed the effect of a matter 
density noise on the MSW effect and found that the presence of noisy matter 
fluctuations weakens the MSW mechanism, thus reducing the resonant conver- 
sion probabilities. These papers, nevertheless, did not take into consideration 
KamLAND data. 

In order to analyze the effect of a noisy density on the neutrino observations, 
we must consider the evolution equations for the neutrino when the density 
is given by a main average proflle perturbed by a random noisy fluctuation. 
This is done starting from the standard Schrodinger equation [35,28]. 

The evolution for the Vf. - Vy system is governed by 
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Here E is the neutrino energy, 6 is the neutrino mixing angle in vacuum, Am^ 
is the neutrino squared mass difference and the matter potential for active- 
active neutrino conversion reads 

Ky{t)^^^p{t){l-Y^), (6) 



where Vey is the potential, Gp is the Fermi constant, p is the matter density, 
rup is the nucleon mass and Yn is the neutron number per nucleon. ^ is the 
fractional perturbation of the matter potential. 

To calculate the survival probability of the neutrinos it is necessary to average 
the terms of the evolution equation over the random density distribution. Here 
we assume a delta-correlated Gaussian noise used in Refs. [35,28]: 

<SAl^+'> =0, < 5A,y{t)5A,y{h) > ^2k{t)S{t-h), (7) 



where the quantity k{t) is defined by 

m = < sA%it) > L{t) = \v!y{t) <e> m, (s) 

where L{t) is the correlation length of the perturbation. The averages, denoted 
by <>, are realized in space for one correlation length, and the delta-correlated 
noise means that the fluctuations in density are completely spatially uncorre- 
cted for separations larger than one correlation length. L{t) must obey the 
following relations 

Ifree « L{t) « A„(t), (9) 



where I free is the mean free path of the electrons in the Sun and A^(t) is the 
characteristic neutrino matter oscillation length. 

In order to guarantee that this condition is satisfied in the whole trajectory of 
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Fig. 1. values of the mixing angle, and several values of the perturbation amplitude, 
C = 0% (sohd line), ^ = 2% (long dashed hue), ^ = 4% (dashed line) and ^ = 8% 
(dotted line). 

the neutrino inside the sun, we assume, following Nunokawa et al. [28], that 

27r 



L{t) = 0.1 Xm{t) = 0.1 



u;{ty 



(10) 



where u(t) is the frequency of the MSW effect, given by 
u;\t)=A{Al{t) + Hr 
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Eqs. (10) and (11) give a value of L{t) between 10 and 100 km in the resonant 
region for the LMA effect, that is, for the inner part of the Sun. 

We calculate the survival probability of the neutrinos solving the evolution 
equation, considering the equilibrium density profile given by the Standard 
Solar Model [26]. The effect of the density noise can be seen in Fig. 1, where 
the survival probability of the neutrinos is plotted as a function of Am? /AE 
for four values of density perturbation, ^ = 0%,2%,4% and 8%. We can 
observe that the effect is amplified when sin^29 = 0.5, rather than when it 
is sin'^26 = 0.9. Considering that the best fit value of Am^ obtained from a 
combined analysis of the solar neutrino data and KamLAND observations is 
of order 10~^ eV^, we can say that the neutrinos with energy around 10 MeV 
will be the most affected by the density noise. 
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Fig. 2. LMA region for different values of the perturbation amplitude, at 95% 
C.L. for several values of the perturbation amplitude, ^ = 0% (solid line), for 
C = 2% (long dashed line), ^ = 4% (dashed line) and ^ = 8% (dotted hue). We 
also present the allowed region for KamLAND spectral data, for the same C.L.. In 
the right-handed side of the figure, the combined analysis of both solar neutrino 
and KamLAND observations is shown. 

In Fig. 2 we present the (sin^ 29, Arri^) parameter space comparing the results 
obtained for solar neutrinos with the allowed regions obtained from Kam- 
LAND observations, for four values of density perturbation, ^ = 0%,2%,4% 
and 8%. We observe that the values of the parameters Am^ and sin^ 29 that 
satisfy both the solar neutrinos and KamLAND observations are shifted in 
the direction of lower values of Am^ and sin^ 29 as the amplitude of the den- 
sity noise increases. In the left-handed side of Fig. 2, the best fit point of 
the solar analysis with no perturbations lies in (sin^ 29, Am^) = (0.82, 7.2 x 
10~^eV^), with a minimum = 65.4. For ^ = 4%, the best fit point goes 
to (sin^ 26*, Am^) = (0.71,3.5 x 10"^eV^), with a minimum = 64.0, while 
for ^ = 8% we have (sin 29, Am^) = (0.33,1.4 x lO'^eV^), with a mmimum 
= 64.2. The numbers of degrees of freedom (d.o.f.) in this analysis is 78, 
obtained from 81 data points, 2 oscillation parameters and ^. 

The right-handed side of Fig. 2 shows the combined analysis involving both 
solar neutrino and KamLAND data. The best fit point of this analysis when no 
perturbations is assumed lies in (sin 29, Am^) = (0.85,7.2 x lO'^eV^), with 
a minimum = 71.1. For ^ = 4%, our best fit point goes to (sin^ 29, Aw?) = 
(0.81, 7.2 X 10~^eV^), with a minimum = 71.5, while for ^ = 8% we have 
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(sin^ 29, Am^) = (0.52, 1.9 x lO^^eV^), with a minimum = 75.5. Here, the 
number of d.o.f. is 91 . 

The main consequence of introducing random perturbation in the solar matter 
density is the appearance of entirely new regions in the (sin^ 29 x Am^) which 
allow simultaneous compatibility of solar neutrino data and KamLAND ob- 
servations. Besides the two standard LMA regions, shown in the right-handed 
side of Fig. 2 by the continuous lines, which were called high and low-LMA in 
Ref. [11], we find a new region displaced toward smaller values of sin^2^ and 
Am^ which we call VERY-low-LMA. There Am^ = (1 ^ 3) x 10~^ eV^ and 
sin^ 29 = (0.4 ^ 0.8), at 95% C.L., obtained for ^ = 8%. 

With no perturbation, the solar data restricts the value of the mixing angle, 
while the KamLAND data can restrict the value of Am^. A combined analysis 

leads to Am^ = (5 ^ 9) x 10"^ eV^ and sin'^29 = (0.67 ^ 0.98) at the low- 
LMA region, and Am^ = (1 ^ 2) x 10"^ eV^ and sin^ 29 = (0.70 ^ 0.95) at 
the high-LMA region, at 95% C.L.. 

Including the effects of random perturbations, the solar data is not able any- 
more to restrict the value of the mixing angle, and both mixing and Am^ 
are more restricted by the KamLAND data. And, since the LMA regions arc 
displaced to lower values of mixing and Am^, the allowed KamLAND region 
with lower values of Am^ becomes important. In fact, such a displacement 
has already been noticed by Nunokawa et al. [28] without considering the now 
available data like as zenith spectrum from solar neutrino observations as well 
as data from KamLAND. Therefore their focus was not only on LMA solution 
but also small mixing angle parameters with active and sterile neutrinos. 

In Fig. 3 it is shown as a function of the perturbation amplitude, minimized 
in Am^ and sin'^29. We can see that even for high values of the perturbation 
amplitude we still can have a viable solution. We notice that even in a noisy 
scenario the compatibihty of solar neutrino and KamLAND results is still 
good. In fact, although the absolute best fit of the analysis lies on the non- 
noise picture where ^ = 0, we observe that (x^ — Xmin) < ^ ^"^^ < C < 8%; 
showing a new scenario of compatibility. 

We conclude the paper arguing that random perturbations of the solar matter 
density will affect the determination of the best fit in the LMA region of the 
MSW mechanism. Taking into account KamLAND results, the allowed regions 
of Am^ and sin^29 at 95% C.L. moves from two distinct regions, the low-LMA 
and high-LMA, when no noise is assume in the corresponding resonant region, 
to an entirely new VERY-low-LMA, when a noise amphtude of 5% < ^ < 8% 
is assumed inside the Sun. 

This represents a challenge for the near future confront of solar neutrino 
data and high-statistic KamLAND observations. If KamLAND will determine 
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Fig. 3. a function of ^, the perturbation amplitude, where we minimize in 

neutrino parameters. Here the number of degrees of freedom is 91. 

sm^29 < 0.5 and Am^ < 4 x 10~^ eV^ it can be necessary to invoke ran- 
dom perturbations in the Sun to recover compatibility with solar neutrino 
observations. 



Note added: After completion of this paper we received two new articles 
which analyzed similar pictures as ours. One is the Rcf. [37] which, differ- 
ently from its first version [31], included KamLAND data. In this paper a 
step-function correlated noise has been considered which allows to relax the 
condition given by Eq. (9). Note that for the LMA parameters, the correla- 
tion length obtained from Eq. (9) does not differ very much from the one used 
in this reference which, therefore, arrive to similar numerical results as ours, 
summarized in Fig. 3. 

Also, a new paper [38] arrived to different conclusions, concerning the quality 
of the fit for solar neutrino data together with KamLAND data in a noisy 
scenario. 
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